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We show how a test of macroscopic realism based on Leggett-Garg inequalities (LGIs) can be 
performed in a macroscopic system. Using a continuous-variable approach, we consider quantum 
non-demolition (QND) measurements applied to atomic ensembles undergoing magnetically-driven 
coherent oscillation. We identify measurement schemes requiring only Gaussian states as inputs and 
giving a significant LGI violation with realistic experimental parameters and imperfections. The 
predicted violation is shown to be due to true quantum effects rather than to a classical invasivity 
of the measurement. Using QND measurements to tighten the “clumsiness loophole” forces the 
stubborn macrorealist to re-create quantum back action in his or her account of measurement. 

PACS numbers: 03.65.Ud, 03.65.Ta, 42.50.Dv, 42.50.Ct, 42.50.Xa 


Making analogies to Bell inequalities pQ, Leggett and 
Garg (LG) [2] proposed a test for quantum behavior of 
macroscopic systems undergoing coherent evolution. The 
resulting Leggett-Garg inequalities (LGIs) aim to distin¬ 
guish a hypothesized philosophical position of macroreal¬ 
ism (MR) from quantum physics, and ultimately to test 
this position against nature. The MR position holds that 
arbitrarily low-disturbance measurements should be pos¬ 
sible, contradicting the Heisenberg uncertainty principle. 

To make manifest the LG ideas, several experiments 
have tested LGIs. Nearly all experiments have used mi¬ 
croscopic systems, including single photons EMU, a single 
photon stored in a macroscopic quantum memory [7], de¬ 
fects in diamonds ElE], nuclear spins eeqhh], and cold 
atoms [14]. See Ref. [15] for a review. To date, two exper¬ 
iments have tested LGIs on macroscopic systems outside 
the single-excitation regime: those by Palacios-Laloy et 
al. [16] and Groen et al. mi- These experiments used su¬ 
perconducting qubits and showed a significant violation 
of a LG-like inequality for weak measurements nam]. 

Because continuous weak measurements record a sys¬ 
tem oscillating between two conjugate variables, they 
perturb both variables during the multi-cycle measure¬ 
ment. This guarantees a disturbance and opens a “clum¬ 
siness loophole” [20]: a macro-realist can interpret the 
LGI violation as caused by imperfect (from the MR per¬ 
spective) measurements. As argued by Wilde and Mizel 
(WM) [20], the clumsiness loophole cannot be closed, but 
one can force the macrorealist to retreat to unlikely sce¬ 
narios in which the clumsiness is imperceptible except 
in the LGI test. WM considered ideal projective mea¬ 
surements, which can be well approximated only for mi¬ 
croscopic physical quantities. Quantum non-demolition 
(QND) measurement is a practical alternative suitable 
for macroscopic quantities. A QND measurement has 


both measurement uncertainty and disturbance (of the 
measured variable) near or below the standard quantum 
limit [2L. Originally proposed to detect mechanical os¬ 
cillations in gravitational wave detectors [22] . a strictly- 
defined QND measurement has been demonstrated in op¬ 
tical [21] and in atomic [23] systems. 

Here we show that, in contrast to previous approaches, 
QND measurements can test a macroscopic system 
against a true LGI, i.e. absent additional assumptions. 
The approach closely resembles the original LG proposal 
and maximally tightens the clumsiness loophole. Us¬ 
ing the collective quantum variable formalism [24] [25] . 
we predict a violation for realistic experimental param¬ 
eters 25. [26] . with the possibility of a straightforward 
extension of the analysis to different macroscopic sys¬ 
tems, up to everyday-life scales E3- This disproves a 
well-known conjecture [28] that LGI violation in a macro¬ 
scopic system requires high computational complexity 
and, thus, that suitable systems are unlikely to exist in 
nature (see also Refs. omi). It also extends previ¬ 
ous studies of QND-measurement-based LG tests [31] [32] 
that found no LGI violation with three-measurement pro¬ 
tocols, and it explicitly shows a fundamental difference 
between temporal (LG) and spatial (Bell) non-classicality 
in the macroscopic limit (33|. Our calculation method 
allows a clear discrimination between incidental distur¬ 
bances from, e.g., spontaneous scattering, and essential 
disturbance due to quantum back-action. The clumsi¬ 
ness loophole can be tightened as per WM, and doing so 
forces the macrorealist to a position strongly resembling 
quantum mechanics. 

Leggett-Garg inequalities —The simplest LGI [2] 
H5] considers two-time correlation functions Cij = 
( Q(U)Q(tj )), for an observable Q measured at times L, 
and giving discrete values Q(t ) = ±1 [2; or else limited 
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FIG. 1. Schematic representation of the Leggett-Garg test. 
An observable Q is measured at different times U (represented 
along the horizontal axis), giving results Qi . Macrorealism as¬ 
sumes non-invasive measurements, with the consequence that 
correlations, e.g., C 13 = (Q 1 Q 3 ), are equal, independently of 
which sequence was performed to obtain them. Si and S3, 
which differ by the presence or absence of Q 2 , give the same 
C 13 in macrorealism but not in quantum mechanics. LGIs can 
detect macrorealism violations using experimentally observed 
correlations. 

to the range Q(t ) G [—1,1] [E] (cf. Fig. [lj. Under MR, 
the correlations obey 

P3 = C21 + C32 + C31 + 1 > 0, (1) 


The atomic system we consider consists of Na spin- 
1 atoms, described by the collective spin vector J with 
components Jj~ = where jW is the total angu¬ 

lar momentum of the Zth atom. The probe light, or 
“meter” consists of pulses of Nl photons described by 
the Stokes vector S^) for the ith pulse with components 
S k = \{ a \i ur) T where a are the Pauli matri¬ 

ces. The system plus meter are described by the vector of 
observables V = (J, ..., S^ n ^), where n is the total 

number of light pulses. 

The initial state is fully x-polarized, i.e. with (J x ) = 
Na and (S x ^) = Nl/ 2. The angular momentum com¬ 
ponents J y ,J z and the Stokes components S y ,S z have 
zero mean and, due to the large numbers, Na rsj 10 6 and 
Nl ~ 10 8 , are Gaussian distributed to a good approxima¬ 
tion. In this sense, we use a macroscopic number of pho¬ 
tons to perform a measurement on a macroscopic number 
of atoms. We can thus describe the full state using the 
average (V) and the covariance matrix [24, [25, [40. [Tlj 

r ij^vVj + vM-Wty). (3) 

Free dynamics under a magnetic field along the x- 
direction produces the evolution 


while a quantum system can show P 3 = — 1/2 in the 
case of a two-level system [ 2 ], and can reach the algebraic 
bound P 3 = — 2 for an infinite number of levels [34, 35] . 
Equation 0 can be generalized to a class of inequalities 
involving n measurements in time J36J, namely, 


Kn 


N (QiQj) + [2 


>0, 


( 2 ) 


where \n\, denotes the integer part of n. Such inequal¬ 
ities correspond to facets of the Leggett-Garg polytope, 
and therefore provide optimal discrimination of non- 
classical correlations [33 38i. 

QND measurement — An ideal QND measurement 
[23 [22] is an indirect measurement of a variable Ps per¬ 
formed by coupling the “system” (S) to a “meter” (M) 
via an interaction Hamiltonian P int = PsPm, where Ps,m 
are the conjugate momenta to Xs,m, respectively. The 
coupling imprints information about Ps on Am, with¬ 
out disturbing Ps, which commutes with P/. A strong, 
direct measurement is then made on Am, providing infor¬ 
mation about Ps- This leaves the system in a state with 
reduced uncertainty in Ps. In contrast the QND inter¬ 
action produces a back action on the conjugate variable 
Xs, increasing its uncertainty. 

Detailed and accurate models of QND measurements, 
including realistic models for measurement-induced dis¬ 
turbance, have been developed for atomic ensembles 
probed by near-resonant light |25j 39]. We use the collec¬ 
tive quantum variables formalism introduced in Ref. [24] 
and expanded in Refs. [25l 140] , The same formalism can 
be applied to other macroscopic systems m- 


(J) ^ <J 9 > = R x {9){ J), 

Tj 1 y r Ja = R x (9)TjR x ( 9) t , 


where 0 = hiBAt is the rotation angle of the atoms in 
the time At given the coupling constant k = — /ip^F be¬ 
tween the atoms and a magnetic field with amplitude P, 
if /ip is the Bohr magneton and g p the Lande factor; Tj 
refers to just the atomic part of the covariance matrix, 
and R x (0) is the matrix describing rotation about the 
J x axis. A measurement consists of passing a pulse of 
light, which is short relative to the Larmor precession 
time, through the atoms. Faraday rotation, produced by 
a “QND” interaction Hamiltonian H[ nt = gS z J z ,imprints 
the instantaneous value of J z on the light, described by 
the relations (in the small angle approximation) 

5(° ut) = 4 in) + gjP Si in) , (5a) 

4 out) = 4 in) + <?4 in) 4 in) , (5b) 


where g is a coupling constant. J z and S z , which com¬ 
mute with P/, are unchanged. 

The required Hamiltonian has been achieved by dy¬ 
namical decoupling [23, 4[2] and by two-color probing m • 
The linear operator relations of Eqs. (5a) and (5b) 
transform the mean and covariance matrix as 


(V) 1 y Mq (V), 

r 1 y MqTMq, 

where the matrix Mq has the form [44] 


Mq 


(A B at \ 


(6) 

( 7 ) 
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FIG. 2. Numerical evaluation of the reduced LGIs as a 
function of 0. (a) K' n for (from top to bottom) n = 3, 5, 7, 9 
in the presence of scattering, (b) Numerical evaluation of 
Kq. The two lower blue (upper green) curves are results with 
(without, i.e. B a t — 0) back action in Eq. 0 - Solid (dashed) 
lines show results with (without) scattering. All plots are 
obtained using the same parameters, taken from Ref. [25] (see 
text). 


The matrix elements B at and B\ represent the back ac¬ 
tion of the interaction on the atoms (J) and light (S), 
respectively, and are practically given by the terms mul¬ 
tiplying g in Eqs. ( |5a| ), ( |5b| ). In a realistic situation, i.e. 
for finite optical depth (OD), we must take into account 
loss and decoherence due to off-resonant scattering of the 
QND probe light. As discussed in Ref. [25], if a fraction 
1 — X of the Na atoms scatters a photon, this alters the 
variances of the quantum components y and 2 as 


rj^x 2 rj + iv A (i-x)(^ + -)iL. 


(8) 


Note that for constant x —1 (i.e., no scattering 
noise) as OD —>• 00 . Repeated application of these trans¬ 
formation rules gives (J) and Tj describing the now- 
correlated system after all pulses have traversed the en¬ 
semble. Through Eq. (5a), the meter variables Sy^ can 
be taken to represent the correlated measurement out¬ 
comes on the system observable J Z (U) at times U. 

Application to LG inequalities —The n outcomes y = 
(i/ij..., y n ) = (Sy 1 ),. • •, Sy n ^) are distributed according 


to the Gaussian probability density function 


Pr(y) = 4; } (y) = 


exp [-(y - - y)] 

yj (27r) n det Ty 


with mean fi = 0, where Ty is the covariance matrix 
describing y [45.. 

As in the original LG article, we generate a dichotomic 
variable Q(U) = sgn(^) [46] . The correlators Cij can be 
evaluated from the 2 x 2 covariance matrix Ty.. = ( ^ ^) 
obtained as the submatrix of Ty describing the measured 
pulses Sy\ Sy\ In particular Cij = (1 — 2a/n) sgn(B), 
where a = arctan(^/AC/B 2 — 1 ). 

Results —Predicted LGI outcomes are shown in Fig. 
[2] where we evaluate sequences with n = 3,5,7, and 
9 measurements. To directly compare these cases, we 
consider a reduced LG parameter K' n = K n / [^J, not¬ 
ing that K 3 = K 3 . We evaluate K' 3 _ 9 , taking into ac¬ 
count decoherence and losses as in Eq.([ 8 |. For simplic¬ 
ity, we consider the case of equally delayed measure¬ 
ments, i.e. a rotation angle 6 between each of the n 
measurements. Realistic parameters are used: g = 10 -7 , 
Na = 10 6 , Nl = 5 x 10 s , and x = exp [—77 Nl], where 
77 = 0.5 x 10 -9 [25]. No violation is seen with an 
n = 3 protocol; a violation is seen with n = 5, but only 
with very low 77, below current experimental values (not 
shown). For n = 7 and n = 9, the LGI violation is achiev¬ 
able with realistic parameters [see Fig. 1 - In Fig.^b) we 
compare the n = 9 case with and without loss and noise 
introduced due to off-resonant scattering. Note that for 
most #, the effect of unwanted scattering is to reduce the 
observed violation. In contrast, for 0 « 7 r, scattering in¬ 
creases the violation, or can create an apparent violation 
that is absent for an ideal measurement (i.e. with 77 = 0 ). 

The above tests involve a large number of correlation 
terms (e.g., 21 for computing K?). We can considerably 
simplify the protocol and reduce the number of measure¬ 
ment sequences by considering just a triple {Q a , Qb, Q c }> 
extracted out of the n- measurement scheme, and the cor¬ 
responding correlators, namely 


N 3 — C a i) + Cb c + C ac + 1 > 0 . ( 10 ) 


We compute the best achievable K 3 optimizing over 
all possible triples and all possible sequences, since from 
a macrorealist perspective possible additional measure¬ 
ments have no effect. The results are plotted in Fig. [3j 
where it can be seen that a violation of Eq. (10) is ob¬ 
tainable, especially around the points 0 = | and 9 = |. 
The optimal sequences of seven measurements for 0 = ^ 
are also depicted in Fig. [3j For an ideal QND measure¬ 
ment we should get C 35 = C 57 = — 1 , with C 37 < 1 due 
to the various discarded measurements made between Q 3 
and Qy that decorrelate the two measurements and give 
rise to the LGI violation. Compared to the seven-point 
measurement of Fig. [2} this protocol shows less violation 
but requires fewer measurement sequences and involves 
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FIG. 3. Three-point LGI violations within longer measure¬ 
ment sequences. Upper and lower curves show K3 versus 6 
for optimized seven- and nine-measurement sequences, respec¬ 
tively. Both plots are obtained for Nl — 5 x 10 8 including 
scattering and loss effects. Spirals illustrate optimal sequences 
in the seven-measurement protocol with delay 0 — for 
which K3 = C35 + C37 + C57 + 1 • Hollow orange circles indi¬ 
cate measurements used to compute correlators in K 3 , filled 
green circles indicate measurements performed but discarded. 


calculation of a simpler correlation function, potentially 
making it more robust in the presence of experimental 
uncertainties. 

Classical versus quantum effects — Within a quantum 
interpretation, we can ask whether the violation of LGIs 
witnesses a genuine quantum effect or whether it is due to 
the classical invasivity (clumsiness) of the measurement. 
There are two ingredients that contribute to the viola¬ 
tion: the scattering and the quantum back action of the 
measurement on J y . The violation around 0 = it can be 
easily explained in terms of the classical invasivity of the 
measurement: measurements at angles /c7r, which should 
be perfectly correlated or anticorrelated, are decorrelated 
due to scattering effects. On the other hand, the quan¬ 
tum back action is a genuinely quantum effect, required 
by the Heisenberg uncertainty principle. 

Our formalism allows us to distinguish between these 
two contributions, and we do so by simulating a QND 
measurement where the effect of the quantum back action 
is “turned off”, i.e. B at = 0 corresponding to J y out ^ = 
J y n ^ in the input-output relations of Eqs. (5a) and (5b). 

The results in this last case are shown in the green 
upper curves of the bottom of Fig. [2jb) . These results 
show that the violation genuinely comes from the quan¬ 
tum back action effect in most of the cases. Scattering 
becomes important only at some specific phases, and is 
responsible for a significant violation only for 0 approach¬ 
ing 7 r. 

Tightening the clumsiness loophole — WM suggest per¬ 
forming, in addition to the LGI test, auxiliary mea¬ 
surement sequences that prove the individual measure¬ 
ments are non-disturbing m , and describe appropriate 


sequences for projective measurements on qubits. Pro¬ 
jective measurements are unrealistic in the macroscopic 
context, however. We now show how even non-projective 
QND measurements can be proven to be non-disturbing; 
see also [23l |44] . 

Consider two identical, non-destructive measurements 
in rapid succession, i.e, with no system evolution in be¬ 
tween. If the statistics of the first and second measure¬ 
ments agree, the first cannot have disturbed the sys¬ 
tem. We illustrate this with linear measurements of J z 
with known gain g , described as S y out ^ = n + gJz^ and 
Jz ° ut ) = Jz in ^ + d , where the random variables n and d 
are the readout noise and the disturbance to J^, respec¬ 
tively. Considering two identical measurements in quick 
succession, it is easily shown that ( S y 2 ^ — Sy 1 ^) = (d) and 
var(5'y 2 ^) — var (Sy 1 ^) = g 2 var(d), which provide ready 
quantifications of (d) and var(d). Both of these approach 
zero for QND measurements in the limit of high OD (the 
“ideal case” of Fig. § m- It is thus possible to demon¬ 
strate to a macrorealist that the QND measurement does 
not disturb J z . 

Combined with the LGI violation, this puts the macro¬ 
realist in a tight spot, requiring some kind of “colluding 
measurements” (in the words of WM) to explain the LGI. 
To remain within the framework of realistic explanation, 
the macrorealist must believe there is a condition of the 
system after the QND measurement. This clearly in¬ 
volves a change in the state but not of J z ; some other, 
orthogonal variable must change. Moreover, this distur¬ 
bance must give rise to a LGI violation, so it must be 
a variable that, in time, rotates into J z . Given that the 
magnetic rotation is about the x axis, J x does not rotate 
into J Z: and the only option is that the measurement 
disturbs J y . Being orthogonal to J^, this disturbance 
does not show up in the quickly repeated measurements 
of the auxiliary sequences, but it becomes visible later 
as the state evolves, leading to the LGI violation. Re¬ 
markably, this macrorealist explanation reproduces pre¬ 
cisely, if qualitatively, the quantum mechanical explana¬ 
tion. The macrorealist must re-invent quantum back ac¬ 
tion to describe the alleged “clumsiness.” 

Conclusions and discussion — We have shown that 
quantum non-demolition measurements allow true 
Leggett-Garg inequalities to be tested in macroscopic sys¬ 
tems. Protocols involving simple state preparation and 
as few as five measurements can violate a generalized 
LGI, and the degree of violation grows with the number 
of measurements. Using covariance matrix simulations, 
we can distinguish violations due to quantum back action 
from violations due to incidental but unavoidable effects 
such as incoherent scattering. We show how QND mea¬ 
surements can be used to tighten the “clumsiness loop¬ 
hole” in a macroscopic system and force the macrorealist 
to a position closely resembling quantum mechanics. The 
LG test strategy described here uses only Gaussian states 
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and Gaussian measurements, and can foreseeably be ap¬ 
plied to very large objects such as gravitational wave in- 

de Groot, T. Picot, G. Johansson, and L. DiCarlo, Phys. 
Rev. Lett. Ill, 090506 (2013) 

terferometer mirrors [27 . 
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